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Abstract On the basis of a dynamical Relativistic Mean Field (RMF) model we study 
the response of a nucleus on the antiproton implanted in its interior. We solve the 
Vlasov equation for the p-nuclear system and show assuming a moderately attractive 
p optical potential that the compressed state is formed on a rather short time scale 
of about 4 10 fm/c. The evolution of the system after p annihilation is simulated 
using the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model. Finally, 
several sensitive observables to the p annihilation in a compressed nuclear configuration 
are proposed, e.g. the nucleon kinetic energy spectra and the total invariant mass 
distributions of produced mesons. 

Keywords p-doped nuclei • nuclear compression • RMF model ■ BUU model 



1 Introduction 

Creation of the compressed nuclear matter in laboratory is a hot topic in the heavy 
ion community for more than three decades. The main interest here is to extract the 
information about the Equation of State of nuclear matter at high densities relevant for 
neutron star interiors. A well established way to produce compressed nuclear matter is 
to collide heavy ions with high energy. However, in this case the produced system will 
be not only compressed, but also quite strongly excited thermally. 

Recently several authors [Il|2][3','4','5' discussed the possibility to compress the nu- 
clear system without much thermal excitation by implanting in it a hadron ( A, K, p, A) 
which has an attractive nuclear potential. It is believed that the antibaryon potentials 
are strongly attractive. The G-parity transformed proton potential is extremely deep: 
Up ~ —700 MeV, while the phenomenological antiproton optical potential is much 
smaller Up ~ -(100^ 350) MeV [HlEll^. The static RMF calculations [SE] predict a 
strong nuclear compression effect by p and A even in the case of realistic potentials. 
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For instance, the central nucleon density of the O core in the ground state p (8) O 
system is (2 4)po, where pg ~ 0.148 fm~^ is the normal nuclear matter density. 

A natural question which arises here is whether the early annihilation of an an- 
tibaryon will kill the effect of compression? A simple estimate of the stopped p life time 
in nuclear matter at the normal density using vacuum annihilation cross section gives 
T ~ 2 fm/c. The phase space volume available in the baryon-antibaryon annihilation 
could be, however, essentially reduced due to smaller effective masses of the annihi- 
lating pair in nuclear matter. This effect can delay the annihilation up to ~ 20 fm/c 

m- 

In this talk we will present our results of dynamical simulations of p-nuclear sys- 
tems to demonstrate the compression dynamics and also to propose observable signals 
of possible p annihilation in a compressed nuclear environment. Sect. [2] contains the 
short description of a model. The numerical results are collected in SectO while the 
conclusions are given in Sect. ID 



2 The model 

Here we give only a very brief review of our dynamical model. The detailed description 
of the GiBUU model applied in calculations and of its RMF extension can be found in 
refs. [5]|9] and in refs. therein. 

In calculations we use the RMF Lagrangian of Ref. 4 , which includes the nucleon 
and antinucleon fields interacting with the isoscalar-scalar (T-meson field and with the 
isoscalar-vector cj-meson field. The NL3 set of parameters was used. This set provides 
a good description of nuclear matter and ground states of finite nuclei. 

The antinucleon couplings with a- and oj-fields were obtained by rescaling the nu- 
cleon coupling constants as g^j^ — ^go-jv and g^jy — —^g^^N, where < ^ < 1 is 
a scaling factor. Exact G-parity transformed antinucleon scalar and vector fields are 
obtained by the choice ^ = 1. The antinucleon potentials consistent with the phe- 
nomenological ones are represented here by choosing ^ = 0.3. 

The dynamics of interacting nucleons and an antiproton was described within the 
transport GiBUU model dealing with a coupled set of semiclassical kinetic equations 
(see [5] and refs. therein). The test particle technique in a parallel ensemble mode has 
been used to solve these equations. 

An antiproton was initialized according to the gaussian density distributions in 
coordinate and momentum space. The centroids of the gaussians were set to describe 
the particle at the center of a nucleus (r = 0) at rest (p = 0). After initialization the 
system was left to evolve according to the kinetic equations without collision terms 
until p annihilation at some preselected time moment tunn- At this time a sudden 
annihilation of the antiproton with the closest nucleon was simulated on the basis of 
the string model JETSET. Once the annihilation is simulated, further evolution of the 
system of nucleons and mesons has been followed by solving the full kinetic equations 
including the collision terms. In this way, we have taken into account, in-particular, 
important processes of pion rescattering and absorption: nN A ^ tvN or tvN — > A, 
AN NN. 
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Fig. 1 Nucleon and antiproton densities (top panels) and potentials (bottom panels) vs 
coordinate z on the axis passing through the center of the p'^'Ca system at selected times 
indicated in the figure. The calculations with the scaling factor ^ = 0.3 (^ = 1) are shown in 
the left (right) panels. 



3 Numerical results 



First, we have studied the compression dynamics of a nucleus under the influence of an 
antiproton implanted in its center, completely neglecting the collision terms in kinetic 
equations, i.e. by taking into account the mean field only. The upper panels of Fig. [T] 
show the axial density distributions for nucleons and an antiproton at different times. 
In the lower panels of Fig. [1] the potentials defined as Uj = gujj(^'^ + gaj<^ {j = N, p) 
are also shown. The nucleon density at the center grows quickly reaching its maximum 
aheady at t = 10 fm/c for the reduced p-couplings and even faster a.t, t — 4 fm/c, 
for the G-parity-motivated couplings. The maximum nucleon density reached in this 
process is p™'''^ — 2po {Pn^^ — 3po) for ? ~ 0.3 = !)■ The compression of the nuclear 
system leads also to the growth of the antiproton density at the center. The system 
stabilizes on the time scale of several tens fm/c. 
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Fig. 2 Time evolution of the nucleoli densities at the centers of the p^^O, p'^^Ca, and p^'^^Pb 
systems for two values of the scaling factor ^ = 0.3 (left panels) and ^ = 1 (right panels). 
The dotted line shows the calculation without p annihilation. The thin solid, thick solid and 
dash-dotted lines show the results with annihilation simulated at various times tann indicated 
in the figure. The dashed lines show the central nuclcon density for the corresponding ground 
state nucleus without an antiproton. 



Compression process is better demonstrated in Fig. [S] where the central nucleon 
density evolution is presented for the three systems: p^^O, p^^Ca, and p'^'^^Pb. For all 
the systems a compressed state is reached within t = 4 10 fm/c and exists up to 
t = 100 fm/c. At later times the stability is gradually lost due to the test particle 
escape beyond the computational grid. 

To understand how compression can possibly influence observables, we have simu- 
lated the annihilation of an antiproton at different preselected times faim- The evolution 
of a central nucleon density for different tann is also shown in Fig. [5] If the annihilation 
takes place in a compressed state, the system expands and, after 30 -f- 40 fm/c, the 
central density drops below pQ. For the light and medium systems, p^^O and p^'^Ca, 
the central nucleon density stays below 0.6po for about 30-^50 fm/c, which is sufficient 
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Fig. 3 Kinetic energy spectra of emitted nucleons in the cm. frame for various pA systems 
and values of the parameter ^ as in Fig. [2] Different histograms correspond to different values 
of the annihilation time tann indicated in the key. 

for the small density perturbations to grow according to the spinodal mechanism [lOj . 
This process should eventually lead to the multiple fragment formation. 

Fig. |3]shows the kinetic energy spectra of nucleons emitted from the system after p- 
annihilation. They are mostly kicked-out by pion rescattering via intermediate Zi(1232) 
resonance: nN — + Z\ -kN [11]. In the case of annihilation in the compressed state, 
the nucleon spectrum acquires a high-energy tail with clearly stronger effect for lighter 
systems. This can be explained by the following mechanism: During the compression 
stage the nucleons in the central zone increase their kinetic energies due to falling down 
into the deep potential well created by the antiproton (see lower panels of Fig. [l|. 
After p-annihilation the potential well disappears suddenly, and the fast nucleons get 
released. This mechanism combined with additional kicks by pions leads to a visible 
enhancement of the energetic nucleon production. 

In Fig.|4]we show the distribution of annihilation events in the total invariant mass 
Minv of emitted mesons. In the absence of any in-medium effects, this distribution 
should be sharply peaked at 2m j^. The final state interactions of mesons, i.e. reseat- 
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Fig. 4 Distributions of annihilation events in the total invariant mass of emitted mesons for 
various pA systems and values of the parameter § as in Fig.E] Different histo grams correspond 
to different values of the annihilation time iann indicated in the key. 



tering and absorption on nucleons, lead to the strong spreading of the distribution 
toward smaller Afjuv Additionally, the mean- field effects reduce the invariant energy 

of annihilating p-nucleon pair as -^5 ~ 2m + Ujy[ + Up < 2m]\f. This further shifts 
the distribution to smaller Afjnv- The effect is especially strong for annihilations in 
compressed configurations. 

The above results have been obtained for the artificial case of an antiproton sitting 
at rest in the center of a nucleus. It is, certainly, more realistic to study the compression 
of the nucleus induced by moving p. As an example, we have considered the central 
p^^C collision at iSiab = 180 MeV. In order to compute the mean field, we have 
performed simulations for 1000 parallel ensembles of test particles. The probability of p- 
annihilation within the time interval At was calculated l-exp{-rp''At), 
where r^"" is the annihilation width and At = 0.2 fm/c is the time step of the 
simulation. In each run, the annihilation time for all parallel ensembles has been chosen 
randomly using Pann- Further details of the model implementation for the p- induced 
reactions on nuclei will be published elsewhere. 
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Fig. 5 Nucleoli and antiproton density and potential profiles along collision axis at different 
time moments for the p(180 MeV)+'^^C collision at the impact parameter b = 0.5 fm. The 
mean field was computed with reduced p couplings (t; = 0.3). 



In Fig.[5]we show the density distributions of nucleons and of the antiproton along 
the beam axis for 20 computed runs. One can see that in 19 runs the antiproton 
annihilation happened at the periphery of the nucleus. And only in a single run the 
antiproton has safely reached the central zone. In this run, the nucleon density has 
been increased up to ~ 2po by moving p. This means that with ~ 5% probability one 
can expect the compression effect even in realistic central p-nucleus collisions. 



4 Conclusions 

To summarise, we have performed the dynamical RMF calculations of the cold com- 
pression of nuclei by an antiproton implanted in their centers. The calculations in a 
pure mean field mode resulted in a fast increase of the central nucleon density up to 
2 3po within the time interval of 4 10 fm/c. Without annihilation, the system sta- 
bilizes in a compressed state, which has been numerically traced until 100 fm/c. The 
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density profiles of compressed p-nuclear systems are close to the ones obtained earlier 
by static RMF calculations [3l|4] . 

Due to in-medium effects, the p annihilation can be suppressed. In this case the 
time needed for a compression can be comparable with the life time of an antiproton 
with respect to the annihilation. Therefore, the compression can, indeed, manifest itself 
in some p-nucleus collision events before the antiproton annihilation. 

It has been demonstrated, that the events with nuclear compression have a proba- 
bility of few percent even in p(180 MeV) + ^^C central collisions. At higher energies, the 
fraction of the compression events is expected to be higher due to smaller p-nucleon 
annihilation cross section. 

In the present calculations we have used quite moderate values for the attractive 
real part of p optical potential. Actual determination of Up requires a very careful com- 
parison with experimental data on elastic and inelastic scattering and absorption of 
antiprotons on nuclei. We propose to perform such measurements at the future FAIR 
facility in Darmstadt. In this respect we find very promising to use the event-by-event 
transverse momentum correlations proposed recently to determine the A optical poten- 
tial [12]. We believe that the same method can also be applied for the determination 
of the p optical potential. 
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